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ABSTRACT 

In many nucleated dwarf elliptical galaxies (dE,Ns), the nucleus is offset by 
a significant fraction of the scale radius with respect to the center of the outer 
isophotes. Using a high-resolution A^-body simulation, we demonstrate that the 
nucleus can be driven off-center by the m = 1 counter-streaming instability, which 
is strong in flattened stellar systems with zero rotation. The model develops 
a nuclear offset of the order of one third of the exponential scale-length. We 
compare our numerical results with photometry and kinematics of FCC046, a 
Fornax cluster dE,N with a nucleus offset by 1'.'2; we find good agreement between 
the model and FCC046. We also discuss mechanisms that may cause counter- 
rotation in dE,Ns and conclude that the destruction of box orbits in an initially 
triaxial galaxy is the most promising. 



Subject headings: galaxies: dwarf — galaxies: evolution — galaxies: individual 
(FCC046) — galaxies: kinematics and dynamics — galaxies: structure — galaxies: 
nuclei 
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1. Introduction 

Dwarf elliptical galaxies (dEs), faint, low-surface brightness galaxies with smooth ellip- 
tical isophotes, are the most common type of galaxy in clusters and groups of galaxies (see 
Binggeli & Ferguson (1994)). In hierarchical models of cosmological structure formation, 
dwarf galaxies form first and subsequently merge to form larger galaxies. Understanding 
their formation and evolution is therefore vital to a complete picture of structure formation. 
There is evidence that dEs constitute a heterogenous class of objects, kinematically (rotating 
versus unrotating, Geha et al. (2003)), chemically (Michielsen et al. 2003) and morphologi- 
cally (Barazza et al. 2002). About half of the dEs harbor a central bright nucleus and are 
called nucleated dEs (dE,Ns). dE,Ns appear to be significantly older than normal ellipticals 
and non- nucleated dEs (Rakos & Schombert 2003); upon tidal stripping, their nuclei have 
been suggested as sources of both the recently discovered ultra compact dwarfs (Phillipps et 
al. 2001) and massive globular clusters like a;Cen (Gnedin et al. 2002). In ~ 20% of Virgo 
dE,Ns, the nucleus is significantly displaced with respect to the center of the outer isophotes, 
typically by ~ 1" (~ 100 pc) (Binggeli et al. 2000). There is a tendency for the displacement 
to increase with decreasing surface brightness but no relation between nuclear displacement 
and any other structural or environmental parameter has been found. 

Various models have been proposed to explain these offset nuclei. In this Letter, we 
investigate whether the lopsided (m = 1) counter-streaming instability can reproduce a 
dE,N with an off-center nucleus. This instability has been known since Zang & Hohl (1978), 
and been studied analytically (Sawamura 1988; Palmer & Papaloizou 1990) and with A^-body 
simulations (Merritt & Stiavelli 1990; Levison et al. 1990; Sellwood & Merritt 1994; Sellwood 
& Valluri 1997). Sellwood & Valluri (1997) do not detect it in systems rounder than E6, 
which is mainly due to their rounder models being stabilized by a higher radial pressure. On 
the other hand, Merritt & StiaveUi (1990) find lopsidedness developing in systems as round 
as El but with negligible radial pressure. Partial rotation only introduces a pattern speed 
in an otherwise purely growing instability; simulations find that the lopsidedness produced 
by the instability is long-lived. 

We use a realistic multi-component A^-body model to generate a lopsided system which 
we then compare with observations. In Section 2, we present the A"-body model and in 
Section 3 we compare it with photometry and kinematics of FCC046, an example of such a 
dE,N with an offset nucleus. In Section 4, we discuss ways in which counter-rotation may 
arise in dE,Ns. 
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2. A/^-body simulation 

We have been using A^-body simulations to explore the counter-streaming instability as 
a mechanism for producing lopsided dE,Ns. We have performed both fully self-consistent 
simulations with live nucleus, disk and halo run on a tree code, and restricted simulations 
with rigid hales run on a grid code; both types of simulations will be presented elsewhere. 
Here we present only one simulation which happens to match FCC046 best, consisting of 
live disk and nucleus components inside a rigid halo. The rigid halo was represented by 
a spherical logarithmic potential with core-radius Tc and asymptotic circular velocity v^. 
The initial disk was modeled by an exponential disk of mass M^, scale-length R^, Gaussian 
thickness and truncated at i?f The nucleus of mass was generated by iteratively 
integrating a distribution function in the global potential (Prendergast & Tomer 1970). 
We used a distribution function of the form f{x,v) oc | [— t?)]"^^^ — [— -Emax]^''^|- Here 
Emax = ^totl^^n); the total potential at the nucleus truncation radius, r^, in the disk mid- 
plane. Initial disk kinematics were set up using the epicyclic approximation to give Toomre 
Q = 1.3. Vertical equilibrium was obtained by integrating the vertical Jeans equation. 
The disk and nucleus were represented by 36.8 x 10^ and 3.2 x 10^ equal-mass particles, 
respectively. After setup, we switched the direction of the velocity for half the particles, 
producing an unrotating disk and a rotating nucleus. We used a quiet start (Sellwood 1983) 
to cancel spurious momenta and angular momenta. 

In our units, where Rd — M{— + M^) — G — 1, which give a unit of time 
{Rl/GMy/'^, we set Zd = 0.2, Rt = 5, rn = 0.78, Tc = 5 and vq = 0.648. The simulation 
was run on a 3-D cylindrical polar grid code (described in Sellwood & Valluri (1997)) with 
Nr X X Nz = QO X 64 X 225. The vertical spacing of the grid planes was Sz = 0.0125. We 
used Fourier terms up to m = 8 in the potential, which was softened with the standard Plum- 
mer kernel, of softening length e = 0.0125. Time integration was performed with a leapfrog 
integrator using a fixed time-step 6t = 0.006. To scale units to FCC046, we note that its 
surface-brightness profile dechnes exponentially with scale-length R^ = 4'.'1 0.4 kpc. We 
estimated the stellar mass of FCC046 at M 1.0 x lO^M©, given a S-band total luminosity 
Lb — 2 X 10^Lb,q and using a S-band mass-to-hght ratio M/Lb ~ 5Mq/Lb,q, typical 
for an old (15 Gyr), rather metal-poor ([Fe/H]= —0.5) stellar population. With this scale- 
length and mass, the unit of time corresponds to (i?^/GM)^/^ f« 4 Myr and the velocity unit 
becomes {GM/Rdfl'^ ^ 100 km/s. 
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Fig. 1. — Panel a: a S-band image of FCC046; panel b: the surface brightness of the A^-body 
model ait — 150. The isophotes are spaced by 1 mag/arcsec^ with the //b = 24 mag/arcsec^ 
isophote drawn in white. Panels c and d: the model's velocity field and the velocity dispersion 
field, respectively. The velocity field perturbations have an m = 2 symmetry and range 
between —0.06 (black) and +0.06 (white), corresponding to ±0.1 x Vdrc- 
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The global lopsidedness of the model, as measured by the amphtude of the m — 1 term 
in a Fourier expansion of the surface density, rose exponentially with time, saturating at 
t ~ 70, for an e-folding time r = 5.67. The angular momentum of the disk and the orbital 
angular momentum of the nucleus increase at the expense of the nucleus's internal (rota- 
tional) angular momentum. By the end of the simulation, the nucleus has lost approximately 
20 % of its initial internal angular momentum. The resulting orbital velocity of the nucleus 
(^nuc ~ 0.01) is, however, only a few percent of the local circular velocity (vcirc ~ 0.35), so 
that dynamical friction against a live dark matter halo is not likely to damp the instability, 
as is also borne out by the live halo simulations. 



3. FCC046: a case study 

FCC046 is an = 15.99 dE4,N in the Fornax cluster. Photometry and long-slit spec- 
tra, with an instrumental broadening cxinstr = 30 km/s, in the wavelength region containing 
the strong absorption hnes of the near- infrared Call triplet, were observed at the VLT with 
F0RS2 in November 2001 (photometry) and 2002 (spectroscopy). Seeing conditions were 
typically O'.'S FWHM. A detailed analysis of the photometry has been presented in De Ri- 
jcke et al. (2003b), while details of the kinematics of FCC046 will be presented elsewhere. 
Additionally, the HST archives contain WFPC2 images of FCC046 in the filters F814W and 
F555W. Availability of all these data was the sole reason for selecting this object. 

FCC046 contains a bright nucleus, which is displaced by 1'.'2 with respect to the outer 
isophotes, or 30% of R^. The lopsidedness involves not just the nucleus, but extends to 
~ 2i?d- Color maps show no traces of large amounts of dust but suggest the presence of a 
blue stellar population. The H« image reveals ionized emission regions, evidence for ongoing 
star formation, within ~ li?d. The position of the blue stellar population does not coincide 
with the lopsidedness of FCC046 and is not likely to be the cause of it (although it may be 
a consequence of it). FCC046 is in the outskirts of the Fornax cluster; its nearest neighbors 
are dwarfs over 100 kpc (in projection) away, and the nearest large galaxy is > 200 kpc 
distant. Interactions are therefore not likely to account for its morphology, unless a weakly 
damped mode has been excited (Weinberg 1994). The kinematics of FCC046 (mean velocity 
and velocity dispersion) were extracted from the spectra following the method discussed in 
Dc Rijcke et al. (2003a). Despite its E4 fiattening, the main body of FCC046 has zero mean 
rotation. The lack of rotation excludes the possibility that the eccentric instability (Adams 
et al. 1989; Miller & Smith 1992; Taga & lye 1998) is responsible for the lopsidedness. 

The nucleus is resolved by HST; it has a FWHM= 0'.'27, significantly larger than the 
FWHM= 0'.'15 of the stars in the same image, and therefore cannot be a foreground star. 
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The nucleus has the same systemic velocity as the galaxy which rules out a chance projection 
of a background object. It comprises ~ 10 % of the S-band luminosity of FCC046 and is 
therefore probably quite massive. We can use the fact that the outer isophotes of the A^-body 
model do not move appreciably and their center coincides with the model's center-of-mass 
(COM) to estimate the relative mass-to-light (M/L) ratios of the disk and nucleus. In order 
for the mean position of the i?-band light distribution to coincide with the COM, the M/ L 
ratio of the nucleus has to be less than half that of the disk. If the disk consists of an 
old, metal-poor population with M/ Lb ~ 5Mq/ Lb,q, this would mean that the nucleus has 
M/Lb < 2Mq/Lb,q, typical for a very young stellar population (age < 2 — 3 Gyr). This is 
consistent with the blue color of the nucleus and its strong Pa absorption (Michielsen et al. 
2003). The rotation axis of the nucleus is not aligned with that of the galaxy's main body, 
showing evidence for minor axis rotation with an amplitude of w 5 km/s. 

The surface density, the velocity field and the velocity dispersion field of the A^-body 
model at = 150 are presented in Fig. 1. The viewing angle and the spatial scale were 
chosen such that the model's surface density gives a fair approximation to the 5-band image 
of FCC046. The nucleus is offset by ~ 1" to the south-west with respect to the outer 
isophotes and, in response, there is an excess in the disk density to the north-east, as in 
FCC046. Clearly, the instability engenders a global response of the whole stellar system. In 
Fig. 2, we compare the surface photometry of FCC046 with the A^-body model (convolved 
with a Gaussian to mimic the seeing conditions of the observations) in greater detail. Simply 
shifting the magnitude zero-point of the A'"-body model and applying the same spatial scale 
as in Fig. 1 sufficed to match the surface brightness profile of FCC046. As in FCC046, the 
isophote centers of the model oscillate along the major axis; in the model this keeps the 
total centcr-of-mass fixed. Just outside the nucleus, the density contours have a distorted 
shape which, in projection, results in the isophotes being flatter (around ~ 5" in FCC046 
and around ~ 3" in the model) than at larger radii. On the whole, the model successfully 
reproduces gross features in the surface photometry of FCC046. 




^ 5 10 15 20 5 10 15 20 

Tp = Vab (arcsec) rp = Vab (arcsec) 



Fig. 2. — Surface photometry of FCC046 (left panels), compared with the A^-body model 
(right panels), a and h are respectively the semi-major and minor axes of the isophotes. In 
the top panels, the i?-band surface brightness (in mag/arcsec^) of FCC046 is compared with 
that of the model. In the next panels, the eUipticity e = 10(1 — 6/a) of the isophotes is 
shown. In the bottom panels, the position of the center of the isophotes is presented. Amaj 
and Amin (in arcseconds) quantify the position of the isophotes along the major and minor 
axes, centered on the nucleus. 
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Fig. 3. — Part of the spectrum of FCC046, around the strongest Call absorption hne. The 
radial bins over which the spectrum has been summed to boost the S/N are indicated in 
the figure (left column, top three panels: along the major axis, towards east; right column, 
top three panels: along the major axis, towards west; bottom panels: along minor axis, 
towards north (left column) and towards south (right column)). Radius is measured from 
the nucleus. Observed spectrum: black line; model spectrum: grey line (sec text for details). 
This absorption line appears split into two distinct lines, separated by 70 ± 25 km/s, in the 
outermost major axis spectra. The S/N in these panels varies from 5-10 to 20-25; adjacent 
spectra from the left and right columns have comparable S/N. The minor- axis spectra have 
the same continuum level as the outermost major-axis spectra. Clearly, the Call line does 
not appear split in the minor-axis spectra. The model reproduces rather well the major-axis 
line width and the observed line-splitting. 
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Fig. 1 shows that, while the model's mean velocity is zero, bi-symmetric perturbations 
of order ~ 0.06 are present. The velocity dispersion profile is asymmetric, rising less rapidly 
towards the left side of the model, in the direction of the light excess. The apparent flattening 
of FCC046 and its zero mean rotation already hint at two equal-mass counter-streaming 
stellar populations. Although at the limit of what can be done, given the spectral resolution 
and the low S/N-ratio of the data, the strongest Call absorption line in the major-axis 
spectrum appears split into two distinct lines (see Fig. 3), separated by 70 ± 25 km/s; 
the other, weaker, lines have much lower S/N-raiios. This absorption line is unaffected 
by the bright sky OH-emission lines. (If a badly removed sky hne were the cause of the 
line-splitting, it would also manifest itself in the center of the galaxy and in the minor axis 
spectrum, which is not the case.) What is particularly compelling about these data is that 
the apparent bimodality is present (and absent) where the model predicts it should be (before 
the signal disappears in the noise), and that we see similar splitting in three different spatial 
bins. To illustrate this, in Fig. 3 we plot the observed spectrum and the model spectrum, 
consisting of a fiat continuum and a very narrow ((5-function) absorption line, broadened with 
the F0RS2 instrumental profile and the model's line-of-sight velocity distribution (LOSVD), 
summed over the same spatial bins as the observed spectrum. We did not use an observed 
stellar template since the width of the individual LOSVD peaks is smaller than one F0RS2 
pixel, making an accTiratc calcTilation of a synthetic galaxy spectrum impossible. Except 
for underestimating the minor-axis velocity dispersion, the model reproduces rather well the 
major-axis line widths and the observed line-splitting. 

However, these observations, while very suggestive, need to be corroborated at higher 
spectral resolution to resolve the two peaks in the LOSVDs. In Fig. 4, we present a mo- 
saic of LOSVDs at different points in the model. Moving out along the major axis, the 
LOSVDs become more double-peaked due to the increasing velocity separation between the 
two counter-streaming stellar populations and the declining velocity dispersion. The small 
velocity-field perturbations of Fig. 1 are reflected here in the LOSVDs being slightly asym- 
metric. Although the details of the LOSVD field in Fig. 4 are model-dependent, an observa- 
tional confirmation of its salient features would lend considerable support to the hypothesis 
that the counter-streaming instability is responsible for off-center nuclei in dE,Ns. 
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Fig. 4. — A mosaic of the LOSVDs at different points in the projected model (the middle 
panel is centered on the center of the outer isophotes). Big and small tickmarks indicate 
a velocity v — ±1 and v — ±0.5, respectively, as indicated below the leftmost panel. The 
isophotes of the model, separated by 1 mag/arcsec^, are plotted in grey-scale. 
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4. Discussion 

If the counter-streaming instability is indeed working in a sizable fraction of the dE,Ns 
to drive nuclei off-center, it is natural to ask why these dEs have counter-rotation in the first 
place. The first detections of counter-rotation in galaxy disks (Rubin et al. 1992; Merrifield 
& Kuijken 1994) lead to the hypothesis that the capture of gas with retrograde rotation 
and subsequent star-formation produced two counter-rotating stellar disks. This hypothesis 
predicts stellar populations of different ages and metaUicities. Generally, the two disks will 
not have identical masses, so some net rotation is expected; this is, therefore, an unlikely 
explanation for FCC046. 

Tremaine & Yu (2000) suggested that counter-rotation is produced when a triaxial 
halo with an initially retrograde pattern speed slowly changes to a prograde pattern speed. 
However, this situation probably does not occur very often. Moreover, this mechanism 
requires evolution on long timescales ~ 10^° yr and hence requires dE,Ns to be old and 
unperturbed. 

A more likely mechanism is the destruction of box-orbits, which are the backbone of 
any slowly rotating triaxial mass- distribution, by a growing nucleus. Stars on box orbits 
come very close to the nucleus and thus can be scattered into loop-orbits; conservation of 
angular momentum then requires that direct and retrograde loop orbits be equally populated 
(Evans & CoUett 1994; Merritt & Quinlan 1998). If the nucleus has a mass larger than 2% 
of the total mass of the galaxy, evolution towards an axisymmetric shape requires a few 
crossing times (~ 10^ — 10^ yr), while our simulations show that the counter-streaming 
instability develops on similar timescales. Late infall of gas (Consehce et al. 2003) (perhaps 
as a result of harassment (Moore et al. 1996; Mayer et al. 2001)) or globular clusters (Lotz 
et al. 2001) may be responsible for growing the nucleus. If this explanation holds, then 
the two counter-rotating stellar populations should have equal scale-lengths and chemical 
properties and nearly equal mass. Additionally, the nucleus may retain some of its initial 
angular momentum; we note that the nucleus of FCC046 is indeed rotating. 

To summarize: using A'^-body simulations, we have shown that the counter-streaming 
instability is a viable explanation for at least the dE.N galaxy FCC046. The models can 
produce offsets as large as observed and successfully reproduce distinctive features in the 
surface photometry of FCC046, such as the oscillation of isophote centers. Although at the 
limit of what can reliably be extracted from the spectra, we find a tantalizing hint of counter- 
streaming in the major-axis spectra of this galaxy, which, however, needs to be corroborated 
at higher spectral resolution. 
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